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I. INTRODUCTION 

Quantum computing offers a new perspective in simulating physical systems [1,2]. While the simulation of a larger 
quantum system is intractable for a classical computer, it should be possible of course on a quantum computer. 
Quantum computers are, from a theoretical point of view, highly simplified quantum systems. A simulation of a 
real physical system on such an arrangement is therefore also interesting necessarily because of the usefulness of the 
simulation result, but also because it may reveal fundamental aspects of the real physical system, which may be 
obscure otherwise. 

Maybe the most general problem of this kind would be the simulation of general quantum operations [3] , the most 
general dynamics a quantum system can undergo. The idea is to represent the initial state of the system on a quantum 
bit array, and supplement it with another array embodying the environment. Subjecting this system to the effect 
of a quantum logic network, and dropping the environment bits, the resulting state of the data should be obtained. 
Naturally, having infinite resources, i.e. arbitrary number of "environment" bits, any operation can be carried out by 
realizing a unitary representation of the process. For a system with a k dimensional Hilbert space, this requires a k 2 
dimensional environment. Lloyd conjectured [4], that a general quantum operation on a system with a k dimensional 
Hilbert-space can be simulated with a fc-dimcnsional environment, which may be in a mixed sate initially. This 
conjecture was falsified by Terhal et al. [5], who give bounds to the dimensionality of the environment required to 
simulate certain quantum operations on a single qubit system. Zalka and Rieffel [6] give a simpler proof via explicit 
construction of counterexamples. More optimal possibilities for the simulation, such as closed loop schemes were 
discussed in detail in Refs. [7,8]. 

Programmable quantum gate arrays or quantum processors [9-12] provide another possible approach to the problem. 
In this case, the parameters of the dynamics should be stored in the initial state of the "environment" (which play a 
role of the "program register" - see below). The quantum processor has two input registers, one storing the state of 
the system to be simulated (the "data register"), while the other contains the data of the operation to be performed. 
After the operation of the circuit, the remains of the program state in the program register may either be dropped 
or be subjected to a measurement. The latter is the probabilistic regime of the quantum processor. The result is 
accepted, if the measurement yields a given result. It was shown by Hillery et al. [11,12], that in this case, any 
operation can be implemented, though some of them with quite low probability of success. Recently Vidal et al. [13] 
have presented a probabilistic scheme implementing unitary transformations with rather high probability. 

In this paper we restrict ourselves to the "deterministic" case, that is, the remains of the program register are 
dropped. Hillery et al. [14] have studied the possibility of implementing general quantum maps in such a scheme. 
They have shown, that there are certain limitations, e.g., an amplitude damping channel cannot be implemented. 
In addition to this, we show that different quantum circuits such as quantum teleportation [15] can be used as a 
programmable quantum circuit for some purpose, as it is related to a phase damping channel [16]. 

We focus here on a specific subset of quantum operations, namely Markovian dynamics. These processes are the 
most important ones in the description of quantum decoherence. The question of simulating Markovian dynamics in 
the quantum computing context was investigated by Bacon et al. [17]. As one of their main results, these authors 
provide a decomposition rule to build more complex dynamics from simpler "primitive" operations. On the other 
hand, Bacon et al. apply the formalism of Gorini, Kossakowski and Sudarshan (GKS) [18], providing an excellent 
framework for the treatment of Markovian processes. We approach the problem of simulating Markovian dynamics 
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differently in a way that is related to programmable quantum gate arrays. We restrict ourselves to single qubit 
dynamics, with a single qubit program register. We intend to simulate the generator of Markovian dynamics on 
a programmable quantum logic array, as the generator characterizes the time evolution completely in this case. A 
single run of the processor simulates an infinitesimal time step. As a further restriction we require that the length 
infinitesimal time step should be encoded into the initial state of the program register. 

The paper is organized as follows: in Section II we review some definitions concerning Markovian dynamics, 
introducing Lindbladians, the generators of the dynamics. The standard GKS matrix-representation introduced 
by Gorini, Kossakowski and Shudarsan in Ref. [18] is reviewed, and a practical recipe is given for finding the GKS 
matrix from of a Lindbladian operator in the one-qubit case, utilizing the affine representation. In Section III we 
describe the general idea of simulating a generator, as understood in this paper. Section IV describes a reversible 
scheme capable of simulating a phase damping channel about an arbitrary axis. A geometrical interpretation of the 
result is also given. In Section V we discuss an application of Bennett's telcportation scheme in this context. In 
section VI the results are summarized and conclusions are drawn. 



II. MARKOVIAN DYNAMICS OF A QUBIT 

In this Section we review the definition of Markovian semigroup and its generators very briefly. The most general 
operation, that a state of a qubit can undergo is described by a completely positive linear, trace preserving map acting 
on the set of density operators, also called a superoperator. This may be written in the Krauss-representation [3] as 

£{Q)=Y J E k QEl (2.1) 

k 

where the E^-s are positive operators such that ^E^Ek = 1. 

In order to introduce Markovian processes, one equips the set of superoperators with a continuous parameter t 
which is called the time. Stationary and Markovian dynamics obey the property 

£ti£*2 — £*i+«2! (2-2) 

with t\,ti > 0. The set of superoperators with property (2.2) form a one-parameter semigroup, the Markovian 
semigroup. We also require the property 

So = 1, (2-3) 

(where 1 stands for the identity superoperator) to hold. 

The property in Eq. (2.2) enables us to define the generators of the semigroup as 

L( Q )= lim tM^l. (2.4) 
w t^o+o t 

The operator L is the infinitesimal generator of the time evolution: 

S t = exp(Li) = lim (i - -l) , (2.5) 



n^oo \ n 

from which follows, that 

dg(t) 



dt 



= L[e(t)} (2.6) 



known as the master equation. 

Let us consider an operator L acting on the Hilbert space of a qubit. The question arises, under what condition 
can this operator represent a generator of the dynamical semigroup. The answer was given by Lindblad [19], and by 
Gorini, Kossakowski and Sudarshan (GKS) [18]. We use the notation of the latter authors. According to this, the 
most general form of a generator of a Markovian semigroup reads 

1 3 

L(g) = -i[H, g] + -^ CijQ&iQ^j] + fa, Qd-j}), (2.7) 
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where the as are the Pauli-matrices. The first contribution on the right hand size describes a possible unitary 
evolution, the Hamiltonian H being a Hcrmitian matrix which can be chosen to be traceless without the loss of 
generality. The second contribution describes the stochastic part of the evolution. The Hermitian positive semidcfinite 
matrix Cjj is called the GKS matrix, and it contains all the information on the nature of the dynamics. Note that in 
order to preserve the trace of the density matrix, Lg should be traceless. In fact, Eq. (2.7) describes the most general 
linear operator taking g to a traceless matrix. 

In the following we describe how to extract the Hamiltonian H and the GKS matrix C, if one is given an arbitrary 
function of a single-qubit density matrix g' = Lg In order to do so we utilize the relation between the GKS matrix 
and the affine representation described also in [17]. 

A generic density matrix can be expanded on the basis of the three Pauli-matrices and the unit operator, obtaining 
its usual real 3-vector representation r displayablc on the unit radius Bloch sphere: 

(2.8) 




so that 



=K i+ fH 



(2.9) 



Direct calculation shows that the real 3-vector L[r] corresponding to the most general Lg in Eq. (2.7) reads 



— h 3 h 2 \ I— 2(C*22 + C33) C12 + C21 C13 + C31 
L\r}=\ h 3 -hAr+\ C 12 + C 21 -2(C n + C33) C 23 + C 32 | r + 2i | C31 - C 1S | , (2.10) 
,-/i 2 fti / V C13 + C31 C 23 + C 32 -2(Cu + C 22 ), 





( C23 


— C*32\ 


r + 2t | 


C31 


-C13 




\c 12 


-C21I 



where 

1 3 

H=-^2hi&i. (2.11) 

»=i 

Thus L appears as an affine linear mapping in the 3 space, the above mapping is called the affine representation of L. 

If we are now given an arbitrary function g', which is traceless and depends on the components of the arbitrary 
density matrix g linearly, we can find the corresponding r[_ according to Eq. (2.8) as a function of the components of rf_ 
representing g. This is a linear inhomogenous operator, which can be always decomposed into a sum of a homogenous 
antisymmetric operator, a homogenous symmetric operator and a vector representing the inhomogenity. For qubits, 
this decomposition of the affine representation of the generator is quite meaningful: according to Eq. (2.10), the 
information on the unitary part of the generator, i.e. the Hamiltonian is encoded into the antisymmetric part of this 
operator, while the real and imaginary parts can be found from the symmetric part and the inhomogenity respectively. 

In the absence of the inhomogenity the generator is zero for the identity operator: 

LI = 0, (2.12) 

thus it describes a unital evolution, which preserves the identity operator. The inhomogenity appears in the complex 
part of the elements of the GKS matrix: for qubits, real GKS matrices correspond to unital dynamics. 

Thus equipped with Eq. (2.10), we have the recipe how to find the standard GKS form of a generator of a dynamical 
semigroup for a quantum bit. 



III. SIMULATION OF INFINITESIMAL GENERATORS 



We intend to simulate the infinitesimal generator Lg of Markovian dynamics, on a single quantum bit. This 
can be understood in several ways. As an elementary step we consider the application of a quantum processor: 
an arrangement of quantum logic gates, and possibly measurements, acting on the qubit in argument, and certain 
ancillary systems. The ancillary systems can be used as a "program register" : their state can influence the action of 
the processor on the qubit in argument, which constitutes a "data bit" in this context. 
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The next question might be, how to interpret the time. A possible generic approach would be to regard the single 
run of the processor as a finite time step At. The repeated application of the processor on the data bit results then 
in a discrete time evolution. One can then examine if this is a stroboscopic version of some valid continuous time 
evolution, and search for the proper master equation. This approach will be discussed elsewhere. Here we adopt a 
simpler interpretation: we expect a single run of processor to simulate an infinitesimal time step: 

Pout = £»in + Lg in dt + 0(dt 2 ), (3.1) 

where dt should be encoded in the l^prog) of the program register. The entire evolution can then be approximated 
with some accuracy by running this process many many times. This implements the equidistant first order Euler 
method [20] of solving Eq. (2.6), but time step, and thereby the accuracy is encoded quantum mechanically. Of 
course, the simulation is completely accurate if dt — > 0, and the number of repetitions tends to infinity. We remark 
here, that simulation of decoherence mechanisms with an array of quantum gates has proven to be fruitful in other 
problems too [21,22]. 

Physically, the simulation scheme can be envisaged as a simple collision model: The data qubit is represented by a 
physical system localized in space. It interacts with flying program bits represented by e.g. spin of particles emerging 
from an oven. Each program bit causes the system to evolve a small time step further. 

It follows from Eq. (3.1), that there must exist a program state, for which dt = 0, and thus g out = Qi n , in the 
quantum processor terminology we would say the processor implements the unit operator. It is a natural requirement 
for this kind of semigroup simulation. Thus our scheme is to some extent similar to the idea of simulating a reservoir 
with beam-splitters of transmittance around unity in quantum optics [23] or a quantum homogenization processes as 
described in Ref. [22]. 



IV. A SCHEME WITHOUT A MEASUREMENT 



First we consider a simple deterministic scheme with a single ancillary qubit, as depicted in Fig. 1. 



qubit 1 
program 



qubit 2 ^ 
input data 




for futher 



processing 



waste 

FIG. 1. The quantum network for the deterministic scheme. A single run of the network simulates an infinitesimal time-step 
in the evolution. 



The first bit contains the program, which is most generally 

|* P rog) = VT^ie 4 x|0) +y/i\l), < e < 1, x real. (4.1) 

We don't consider mixed program states here, we want to study how decoherence associated with Markovian processes 
appears purely quantum-mechanically, and a mixed program state would imply some prescribed classical stochasticity. 
The input data is in bit 2, in the state Q[ n . After the operation of the processor, the contents of the program register 
are dropped, and bit 2 contains the output 

ft>ut - Tn (t>(|* prog )(* prog | ® g in )U^ (4.2) 

which is passed for further processing. 

To have the identity operator implemented, there should be a program state for which g out = g- m holds. We chose 
this to be the program state |0). The most general processor that is possible under such circumstances is a controlled 
U gate: 
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The lower right corner is an SU (2) block: 



t cos(l)e-^ -sin(|)e-^\ 

^ " ^ sin(f y** cos(f)e^ J' (44) 

where we use the standard Euler angle parametrization in the y convention [24] . 

Now we can evaluate Eq. (3.1) with a generic input data state of Eq. (2.9), and program of Eq. (4.1), to obtain the 
effect of a single run of the processor. Note that due to the orthogonality of the program states corresponding to 1 
and U~2 (which is a consequence of the unitarity of U), Eq. (4.2) simplifies to 



Bout 



{l-e)g in + eU 2 QinUl (4.5) 



Thus from the point of view of the effect on the input state, the quantum circuit does nothing else than apply the 
unitary operation U2 on the input state, with the probability e. The difference is that if the operation is realized by 
the quantum circuit, then the information which disappears from the system will be stored in the dropped quantum 
state of the program register. In the other case the information will be entirely classical, one bit per infinitesimal 
timestep: we can be aware, whether the operation U2 was carried out or not. Note that in both of the cases the 
process is reversible, provided that we have the appropriate information at hand. 
According to Eq. (4.5), we can write 

Qout = Qin + 02QinUl ~ Sin)S- (4.6) 

Thus we identify e = dt, and comparing with Eq. (3.1) we get 

Lg in = U 2 ghM - Qin- (4.7) 
The transformation on the Bloch ball corresponding to U2QU2 1S a rotation of the vector r corresponding to g 

r_' =ll(e,<j>,i>)r, (4.8) 
where 1Z is the appropriate element of the adjoint representation of SU(2): 

(— sin cp sin ip + cos 9 cos cp cos ip — cos 9 cos cp sin ip — sin cp cos ip sin 9 cos cp\ 
cos 9 sin cp cos ip + cos cp sin tp cos (p cos ip — cos 9 sin <p sin tp sin#sin</> £ SO(3) . (4.9) 
— sin 9 cos ip sin 9 sin ip cos 9 J 

For the generator in Eq. (4.7) we have thus 

L[r] = (ll(6, ( p,iP)-i)r. (4.10) 

This can be compared with Eq. (2.10) to extract the properties of the generator. The transformation in Eq. (4.10) is 
homogenous, thus the generated dynamics is unital. 

As for the coherent part of the evolution, we get the Hamiltonian 

2£ _l«\ n f)( ( ,-i4> 



sin(0 + ip) cos 2 § — |sin6'( 



V| sin 9 (e** + e-'*) - sin(0 + 0) cos 2 § ) [ ' 

This is zero if cp + ip — (2k + 1)tt or 9 = tt + 2kir holds, which is equivalent to trC/2 = 0. Thus in case of traceless 
unitaries, we obtain a purely stochastic evolution in the sense that it lacks the Hamiltonian part. 
Evaluating the symmetric part of Eq. (4.9), and comparing with (2.10), we obtain the GKS matrix 

sin 2 I sin 2 i sin 2 | sin(-0 — cp) \ sin 6>(cos cp — cos -0) N 

Ce,^=\ \ sin 2 § sin(V> - cp) sin 2 § cos 2 ^ ^ sin0(sin ^ + sin^) I (4.12) 



2 om 2 <""VV Y) 2 2 4 

i sin #(cos cp — cos ip) \ sin 6>(sin cp + sin ip) 



2 6 -2 <p+ip 



cos % sin 



2 ^iii 2 



Specifically, if U2 is traceless because ip = tt — cp holds, we obtain from Eq. (4.12) the following GKS matrix: 

sin 2 I cos 2 cp i sin 2 | sm(2cp) i sin 9 cos < 
Ce,<t>= \ \ sin 2 § sin 2cp sin 2 | sin 2 cp isin#sin</>] (4-13) 
\ sin 9 cos ^ sin 9 sin cos 2 f 
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The matrix in Eq. (4.13) can be interpreted as follows. Consider a unitary operator U £ SU(2) acting on the qubit's 
Hilbert space, and a superoperator £, which describes unitary evolution described by the GKS matrix C. According 
to Ref. [17] unitary conjugation of a superoperator £, that is, 

£' = U^£U, (4.14) 

where U(g) = UgU^ yields another superoperator describing Markovian dynamics as well. The resulting GKS matrix 
is 

C' = llC1l T , (4.15) 

where TZ is the element of SO(3), the adjoint representation of SU(2) corresponding to U, and the T stands for 
transposition. Thus TZ is a real 3-rotation, which can be visualized as a rotation in the Bloch-sphere picture. The 
effect of unitary conjugation is to apply the same operation on a transformed basis. In the actual case, the matrix in 
Eq. (4.13) can be rewritten as 



C M = ft(0/2,^) ft T (0/2,&VO- (4-16) 




The matrix between the two rotations defines a phase damping about the z axis. Thus the process is a phase damping 
about an axis pointing towards the spherical polar angles 0/2, <j). In fact the rotation {J 2 moves the z axis to the 
direction described by the spherical polar angles 6/2, </>, so the direction of the phase damping is "half way" between 
the z axis, and its transform by U%. Note that the angle ip is irrelevant in Eq. (4.16), as it does not influence the polar 
angle of the rotated z axis. The phase damping channel has a single direction which is special, thus it is necessarily 
described by two parameters. In the case if U2 is traceless because 9 = n + 2kir is satisfied, we obtain a phase damping 
channel about an axis in the xy plane. We can conclude that the controlled-U gate with a traceless U2 is capable of 
simulating a generator of an arbitrary phase damping channel in our scheme. 

Returning to the generic GKS matrix in Eq. (4.12) we find that rankCg^^ = 1, thus this general GKS matrix also 
describes a phase damping channel, physically the same process as in the traceless case. The only difference is that 
the evolution is now accompanied by a coherent part, generated by the Hamiltonian in Eq. (4.11). 

V. CONTROL VIA TELEPORTATION 

In this Section we briefly describe another simulation scheme based on Bennett's quantum teleportation. It is 
depicted in Fig. 2 



Bell state detector 

i 



p+L[p]dt 



FIG. 2. Bennett's teleportation scheme as a programmable circuit for Markovian decoherence 

The initial state g impinges at the input of a teleportation arrangement. Two additional qubits are prepared in an 
entangled state 



|* P rog) = Vl - ~e\B ) +Vi(ai|Bi) +a 2 \B 2 ) + a s \B 3 }) , (5.1) 
with I ai 1 2 + I CK21 2 + I CK3 1 2 — 1, which serves as the program state. We use the notation 
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|£ ) = -L(|0>|1>-|1>|0» ; 
|B 1 ) = -L(|0)|0}-|1)|1)) ; 
|B 2 ) = -L(|0)|0} + |1)|1)) ; 

|B 3 ) = -^(|0)|1) + |1)|0)). (5.2) 

for the elements of the Bell basis. Then the usual Bennett teleportation is carried out: a Bell state measurement on 
the two appropriate qubits is carried out, and depending on the result, the appropriate unitary transformation <7j (the 
identity operator or one of the Pauli operators) is carried out. 

For e = we have l^/prog) = I -Bo), the state g is simply teleported: the identity operator is implemented. However, 
for nonzero e we obtain 

g' = (1 — e)g + e (\ai\ 2 aigai + \a.i\ oigoi + \oe3\ 2 a 3 ga 3 ) (5.3) 

as a "teleported" state. Thus the scheme is essentially equivalent to a random application of the Pauli operators. 
Setting dt = e as in the previous Section, we can define 

Lg = \ai\ 2 &igai + |a 2 | 2 <T2£ , 0'2 + |«3| 2 o'3£ | o'3 - Q- (5-4) 
A straightforward calculation shows that the corresponding GKS matrix reads 

Kl 2 \ 

C(ai,a 2 ,a 3 ) = ( |« 2 | 2 . (5.5) 

My 

We have obtained a GKS matrix of rank 3, describing a generic Pauli channel. In this case we have the parameters 
of the dynamics encoded into the program state, too. 

It is worth noting here that this second scheme is indeed irreversible. However, the same process could be simulated 
in a similar reversible framework as in Section IV, utilizing two ancillary qubits, and the universal programmable 
quantum gate array in Rcf. [11,12]. 



VI. SUMMARY 



We have investigated quantum computational schemes for the simulation of a generator of Markovian dynamics on a 
qubit, where the value of the infinitesimal time step is encoded in a quantum state at the input of the device. We have 
described the capabilities of the arrangement if the ancilla to be applied is restricted to a single qubit by characterizing 
the most general Lindbladian that can be simulated. We have found that under these restrictions, a phase damping 
about an arbitrary axis can be simulated. We have also investigated a the Bennett quantum teleportation scheme as 
a possible quantum circuit for performing a similar task. 

The inevitable advantage of this method compared with a classical simulation of quantum dynamics is that it works 
for any, even unknown initial state, which may emerge as an output of another quantum calculation. 

The analysis can be carried out for systems ancillas of larger size but in that case, the number of parameters grows, 
and the problem is less transparent. We remark that for instance, the two-qubit program state enables non-unital 
dynamics as well. 

We believe that the study of simple quantum systems such as those described here facilitates the understanding 
decoherence in general. 
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